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Abstract 

The Hamilton-Jacobi formalism for fermionic systems is studied. We derive the HJ equations 
from the canonical transformation procedure, taking into account the second class constraints 
typical of these systems. It is shown that these constraints ensure the consistency of the solution, 
according to the characteristics of fermionic systems. The explicit solutions for simple examples are 
computed. Some aspects related to canonical transformations and to quantization are discussed. 
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Introduction 



The basic property of fermions, half-integer spin, is a property of microscopic particles 
and is described by Quantum Mechanics. However the classical mechanics of fermions has 
been a very useful tool in their study, in particular regarding their path integral formula- 
tion m. Their Lagrangian and Hamiltonian formulations are well known, in particular the 
Hamiltonian requires of the Dirac formalism of singular systems yj, as far as they have 
second class constraints. The understanding of the Hamilton Jacobi (HJ) formalism can 
be helpful to understand better quantum systems. It has been applied to singular systems, 
as is the case of general relativity in the case of the Wheeler-DeWitt equations. In this 
case the HJ formulation amounts to set the constraints of the Hamiltonian formalism as 
differential equations on the wave function j^. For the WKB approximation of singular 
systems, in particular fermionic ones j^, it would be also useful to have a systematic way 
to obtain the HJ formulation. Here an effort is made with the aim to understand this issue 
better. Although fermionic systems have been widely studied, their HJ formulation has 
been less studied. It has been worked out for bosonic constrained systems jgl, where the 
constraints, obtained from the Hamiltonian formulation, are written as separated equations. 
For fermionic systems, a technique to handle a particular problem in the HJ formalism was 
proposed in [5i]. However, the most general situation is not discussed. In ^ the Giiler 
formalism jol is generalized to include fermionic variables, by the application of the usual 
concepts of classical analysis, whose properties, nevertheless, are not of general validity when 
dealing with nilpotent quantities. In this formulation all the constraints of the Hamiltonian 
formulation are kept as additional equations to the actual HJ equation. In ^ a formulation 
is given, which strongly relies on the Hamiltonian formulation. 

Here we give a formulation of the HJ equation for fermionic systems, which is obtained 
as usual for bosonic theories, from the variation of the action in canonical coordinates, 
co.iden„. the t.a.fo.n,at.n to constant new coordinate. 0. A. . the GUle. to..aHs.. 
we apply the "second class" constraints characteristic of fermionic theories, as additional 
equations. We show that these equations have two important consistency consequences. 
First, from the way we obtain the HJ equation, there are two integration constants for each 
fermionic degree of freedom, and we get a set of equations among these constants, which 
reduce their number to half, as it must be for a first order theory. Further, related to this 
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last fact, as noted in [10|, boundary conditions have to be added to fermionic actions. This 
means that also the generator functions of canonical transformations must satisfy boundary 
conditions. It is shown that the mentioned equations ensure that these boundary conditions 
are automatically satisfied. 

In order to verify the validity of the resulting equations, we consider two examples of 
simple fermionic systems. The solutions to these equations are found to be the same as the 
solutions of the Euler-Lagrange equations. 

In the first section the Lagrangian and Hamiltonian formalisms for bosonic and fermionic 
systems are reviewed. In the second section the HJ equation for fermionic systems is given. 
In the next two sections illustrative fermionic systems of one and two variables are considered. 
In the framework of the second example, in the next section, the relation to the canonical 
transformations is established in a more precise way and, in the last section, the relation to 
the quantum theory is discussed. 

Fermionic Mechanics 

Let us consider a classical system described by n bosonic, even Grassmann, degrees 
of freedom q = {qi,q2, ...,qn) and fermionic, odd Grassmann, degrees of freedom ip = 
{ipi,4'2, ■■■I'^n)- These variables obey the relations 

mj - QjQi = i,j = 1,2, , 



-ipa-ipf^ + -ipf^-ipa = = 1,2, . 

In this case, the Lagrangian function depends on the g's, on the ip^s, and on their respective 
time derivatives 







(1) 



L = L{q, q, ip, t) = L{Q, Q, t), 



(2) 



where Q = (g, 

If we variate the corresponding action 




(3) 



then, in order to get the Euler-Lagrange equations. 





(4) 
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suitable boundary conditions must be imposed. For the bosonic variables it can be done as 
usual by fixing each of them at both extrema. However, for each of the fermionic degrees 
of freedom, as far as they are first order in the velocities, only one boundary condition can 
)e fixed. In this case, for consistency, suitable boundary terms must be added to the action 
Q. For example, the fe^.ouic Uinefc te™ is L,.„ ^ the,, the corrected 

action is given by |10|], 

^-^^/"^^a(tl)^/3(t2), (5) 

with the boundary conditions S[ipa{ti) + i'a(t2)] = 0, that is ^ipaih) = —'^Paiti) + ^a, where 
are constant anticommuting quantities. 
The Hamiltonian is given by 

H{Q,P,t) = qp + i:7i-L = QP-L (6) 

were pi = dL/dq\ vr" = dL/dipa and P = {p, vr). 

The Lagrangian Q is first order in the fermionic variables, i.e. the kinetic term is linear 
in fermionic velocities and the potential does not depend on them. Therefore there are 
primary constraints, 

0« = vr°-r(g,V^), (7) 

where /"(q', '?/') are odd Grassmann functions. In the Dirac formalism for constrained sys- 
tems, these constraints turn out to be second class. Further we suppose that there are no 
more constraints. Thus, due to the fact that in the Hamiltonian © the term ^/'q.tt" com- 
pensates the corresponding kinetic term in the Lagrangian, the canonical Hamiltonian does 
not depend on the fermionic momenta, 

H = H{Q,p,t), (8) 

Therefore, if the Lagrangian is purely fermionic, the Hamiltonian will be given by the po- 
tential. 



HJ formalism for Grassmann variables 



In order to find the HJ equation for the preceding system, we consider the variation of 
the action 

dt + BT, (9) 



S= L{Q,Q,t)dt + BT = r \QP~H{Q,p,t) 
Jti Jti L 



where the second class constraints ((Tj) and, as mentioned in the preceding section, suitable 
fermionic boundary conditions {BT) added to the action, insure us that the variation of 
the right hand side gives the correct equations of motion. For instance, if we consider the 
action L = ^g"'^ipa'4^p ~ ^(^); then we have vr" = —^g°'^'4'p ^ii^l if we impose the boundary 
conditions +ipa{t2)] = 0, then. 



6S = 6 



V„7r° - H{^,t)] dt + '-Mti)'^'"{t2] 



t2 



dH 



*2 / • dV 



0. 



(10) 



Thus, the physical phase space is {2n + /i)-dimensional hyperplane V, solution of ((Tj). 
Let us do a canonical transformation of coordinates, {Q,P) {Q,P). In this case the 
constraints ((Zj) will transform to some constraints 



(11) 



To obtain the HJ equation a variation of these actions is done. 



6S = 6 



6S' = S 



'QP-H{Q,p,t) 



dt + BT 



t2 



QP-H'{Q,P,t) dt+{BTy\ = 



(12) 
(13) 



The relation between integrands is. 



dF 



QP - H{Q, p, t) = QP~ H'iQ, P,t) + — + K, 

dt 



(14) 



where K = and F is a function, whose dependence on the phase space coordinates 

and on time must be such that its variation at the boundary satisfies 6[F{t2) — F{ti)] = 0. 
If now F = F{Q, P, t) - QP, then 



dF -dF ~dF dF d 



hence 



(15) 



~ ~ ~ ■ dF ~dF dF d ~ ~ 

QP - H{Q,p, t) = QP- H'{Q,p, t) + Q— + P^ + — - -{QP) + K. 



A factorization of this gives 



» dF\ dF 



H' + K 



(16) 



(17) 



and additionally the constraints ((Tj) and pij). The sign in the middle term corresponds to 

the interchange of Q and P. Even with these constraints, the quantities Q and P can be 

taken as independent, and we get, 

dF ~ dF dF 

P=_ Q = (-ir-^. H' = H^-^K. (18) 

If, as usual, the new coordinates, P = (p, vf ) and Q = {q,i'), are assumed to be constant, 
which is guaranteed if H' = or, what is the same, H' = K, then the HJ equation will be 
in fact given by a system of equations. If the first equation in fITHj) is applied to the last one 
and to ((Zj), we get the HJ equation. 



H 

as well as, 



OF 

Q,—{Q,P,t),t 



+ |^(Q,P,t) = 0, (19) 



^^^^'^'^l^r(Q). (20) 



Additionally we have the second equation in (jl8p . which can be written as 

dF{Q,P,t) 



dpi 
dF{Q,P,t) 



= even Grassmann constant, (21) 
—tpa = odd Grassmann constant, (22) 



plus the /i (unknown) constraints (jllj) . which eliminate half of the fermionic constants vr). 
Usually, the configuration space variables can be obtained from equations (PT|) and (j221); as 
functions of two integration constants, and this will be the case of ()21|) . from which the 
bosonic variables q can be obtained in terms of p, q, ip and tx. However, before solving the 
equations ((221), we can solve the equations (j^ . Indeed, the fact that (fTT|) are second class 
means that /° are invertible, and a solution for ip in terms of p and vf can be obtained, after 
substituting q by its solution. If this solution is then substituted in (j^ . yU relations among 
the constants vr and ip arise, which will eliminate half of them. 

Consistently with these results, we have that, for an action with standard kinetic fermionic 
term, as a consequence of (j20|) the boundary condition for F will be fulfilled, 

dF 

SF{h) = 5^«(ti)-— (ti) = SMti)nti) = SMt2)nt2) = SF{t2). (23) 

Thus, all these equations (|19H22p . must be solved to get the complete solution for the 
Hamilton principal function (Hpf), 

F{Q,P,t) = S{Q,P,t) + a, 



as well as the solution for the configuration space variables (g, V'); depending on the correct 
number of integration constants, two for each bosonic degree of freedom, and one for each 
fermionic degree of freedom. 



A system L = tpip 

In this section, a simple instance is solved to show the problems which appear when 
fermionic variables are present. 

Consider a system characterized by the Lagrangian L = ipip, the Euler-Lagrange equation 
is = 0. The canonical momentum to the fermionic variable ip, is given by tt = dL/dip = 
—ip. The Hamiltonian is given by Hq = ipir — L = ipn — ipip = ip{'K + %p), where the velocity 
ip can be handled as a new parameter. It vanishes, weakly, according to the second class 
constraint. 

Now the HJ formalism is applied, by substituting vr = dS/dip. In this case the Hamilto- 
nian vanishes and the HJ equation is given by 

where the action depends on the configuration variable ip and on one constant fermionic 
parameter p, i.e. S = S{ip,p). We have as well the equations, 

dS , dS ^ 

where /3 is a constant grassman parameter. 

Due to the fact that the action is bosonic, it must have the form 

S = a{t)p^. (26) 

Applying to it the first equation in (P^j). we get ip = ap, then we apply the second equation 
and (3 = a%l) = a? p. Thus a is a constant, as would result also from (j24j) . Thus, the constant 
fermion solution of the Euler-Lagrange equations turns out. 

Interacting system 

In this section, an interacting system with two fermionic variables ipi and ?/'2; such that 
each one are the complex conjugated from the other ipi* = ip2, will be discussed, 

L = i{ipii)2 + ^2^i) + (27) 



The Euler-Lagrange equations are given by, 



#1 + = 0, - = 0, (28) 



with solutions 



Mt) = ^le^^'/')*, Mt) = 6e(-^'/')*. (29) 
The Hamiltonian of this system is given by, 

H = -ktlJiij2, (30) 

which must be accompanied by the second class constraints, as definitions of the momenta, 
TTi = —itp2 and 112 = —iipi- As a consequence, the Hpf will be the solution of the following 
system of equations, 

HW + '-^^^O, (31) 

z^2 + = 0, (32) 

#1 + 77^ = 0, (33) 

A - ^ = 0, (34) 
dpi 

P2-^ =0, (35) 

Op2 

where, pi = tTj and = ipi are constant odd Grassmann quantities, which satisfy p* = —p2, 
Pi = 1^2) and ipiT^i + ip2'^2 is real. Seemingly, there are too many constants for a first order 
system. However, as will be shown further, the role of the equations is precisely to 

establish relations, corresponding to the second class constraints, between them. 

In order to solve this system, we write the most general even Grassmann function of the 
odd Grassmann quantities pi, P25 V"!? '4^2- 

S{ip, p, t) = So{p, t) + S'i(p, t)V^i + 5'2(p, t)V^2 + S^{p, t)V'iV'2, (36) 

where the fermionic functions are given by Si{p,t) = si(t)pi and S2{p,t) = S2(t)p2, and the 
bosonic ones by So{p,t) = so(t) + soi(t)piP2 and S'3(p,t) = S3o{t) + Ss{t)pip2. 

From the reality of S, we get that the coefficients Sq, Sqi, S30 and S3 must be real and 
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Further we have the conditions 

dS 



-Sipi + (S30 + S3pip2)lp2 = -^^2, (37) 



dS 

TT— = -S2P2 - (S30 + S3PiP2)^l = -^^1, (38) 

dip2 
OS 

TT- = S01P2 + Sii>i + S3p2^1^2 = A, (39) 

dpi 
dS 

TT- = -SOlPl + S2i^2 - S3PlV'lV'2 = ^2- (40) 
Op2 

The first two equations, can be solved for ipi and 4'2, thus obtaining, 

V'l = ^P2, i^2 = . pi, (41) 

which substituted in the second two, (jH^ HOj) . give us, 

A = f^oi - -^^) P2, (42) 

P2 = (-soi + ^^)pi. (43) 

These equations could be identified with the second class constraints, letting only two free 
constants, as corresponds to fermionic theories. 

Taking into account the fact that P2, Pi and p2 are constant, we see that the coefficients 
in (H^ and (PHj) are themselves constant, that is 

soi : = A, soi — = A , (44) 

S30 S3o + t 

if we set A = M — iv, we get now 

soi = VS30 + U, (45) 
S1S2 = v{slo + l). (46) 

Now, in order to write the HJ equation, we note that it has to be written before substi- 
tuting (P^j) and ipUl) into the Hpf, because the time derivative in the HJ equation does not 
act on the fermionic variables ip. Thus we have, 
dS 

+ H = so + S01P1P2 + sipi^i + S2P2V^2 + iho + S3pip2)ipii'2 - ki)ii)2 = 0. (47) 

at 



Taking into account PT|l we get. 



So + 2 \ 1 ■^Ol(s3o + 1) - SiS2(s3o + i) - SiS2(s3o - i) + SiS2(s3o + k) pip2 = 0, (48) 
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or, equivalently sq = 0, and 

(Sgo + l)Soi - S3o(SiS2 + S2S1) + i{SiS2 - S2S1) + SiS2(s3o + fc) = 0. (49) 

From the fact that Si is the complex conjugated of and S30 is real, and writing in (|46p 

V = aa*, we get 

si = a*{s3o + i)e'^, S2 = a(s3o - z)e"*^. (50) 
These equations, together with (pSj). substituted back into (jl^ . give 2f + k = 0, i.e. 

r = -^t + c. (51) 

Thus, if we set ^ = —ae~^'^p2-, we obtain 

V^i = ee^*^* (52) 
^2 = Te-^'*, (53) 



which coincide with the solutions ()29|). 
Therefore, the Hpf is given by 



•S* = So + soi - — P1P2 = So -V'iV'2, (54) 

V sin — I aa* 



where Sq is constant. 

Note that the undetermined functions S30 and S3, do not appear neither in the solutions 
(1521 nor in the Hpf. This can be understood from the form of the Hpf fl36|) and the 
equations (HH), as the terms containing these functions vanish identically. 

Canonical transformation point of view 

Let us consider the solution to the equations (jH^ HOj) for ipi and iIj2- We can get first ipi 
from (jnni); then we substitute it in !^40\i. from which we get, 

i^i = \ (siPi - S01S1P2 - —92(^1(^2 - ^^^PiP2/5i) , (55) 

Sf V S2 S2 / 

^2 = 4 f S2/32 + S01S2P1 + — Pl/?l/?2 - ^^PlP2/?2) • (56) 
si \ Si Si / 

These equations can be written in the form of canonical transformations. In order to see 
it, taking into account the last observation of the preceding section, let us set 6*0 = 0, and 
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S30 = 0, in this case ()55I56|) are given by, 

i>i = — P-P2P1P2, (57) 

si sis2 

^2 = - + (58) 

Note that these equations are symmetrical under the interchange ipi ^ "^2, Pi ^ P2 and 
Pi ^ 132. If we define ip^ = s^Vi, ^2 = S2^p2, t^" = Sipi, vr^ = -S2P2, a = (siS2)~^S3, 



and the function G = 7r°7r2V'i^25 equations ()57p58|l can be written as 



dG 
dG 

^2 = ^2" + «^. (60) 



Similarly, the momenta ()37I38|) . can be rewritten as 

dG 

dG 

In vectorial notation, these equations can be expressed by a single equation 

dG 

Au = u - u° = aJ— , (63) 
du 

where u = (V^i, ^/'2, 712), u° = (V^i, ^^25 ^i; '""2); J is the corresponding Jacobian 
matrix. 

As it can be observed, the function G plays the role of the generating function of a 
canonical transformation (j63j) . with a finite parameter a. The solution needs additional 
conditions, for example "initial conditions" tp" oc vTg and 1^2 ^i- This way to write the 
solution to the HJ equation, could be useful for computing the Van Vleck determinant for 
supersymmetric theories 

Quantum Mechanical features 

Defining, = {ipi,ilj2), "ip" = {4'i,ip2)y ^^e equations (jSZlEHl) can be rewritten as follows 

ij = ip°exp{aS). (64) 
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Indeed, taking 5*0 = as in the preceding section, 



hence 



V'l exp (aS) = + a7r°Vi^2 = 



(65) 



Further, considering that 



~dt 



Sl S2 



(66) 



it can be seen that ip is a solution of the following first order partial differential equation 



1 dijj 
a dt 



a 



J: + -S-H]iP 



a 



where 



E = - 

a 



Sl 

S2 







(67) 



(68) 



S2 



J 



A particular case turns out when the permutational symmetry in eqs. (jFTl is broken by 
the application of second class constraints (jUJ. If moreover eqs. fl^Uj) are applied, we get 



Sl 



k 



S2 
S2 



In this case S can be written as, 



1 

a 



( k 
I— 
2 









k 

'2) 



2a 



1 
-1 



2a 



(69) 



(70) 



If for simplicity S3 is assumed to be a constant, eq. ()67|) is rewritten as 

a dt \ 2a I 



(71) 



This equation resembles the Schrodinger equation. Note that, due to the nilpotency of 
fermionic degrees of freedom. Hip = 0. However, if the Lagrangian (j27|) would be extended 
to a supersymmetric one, by the addition of two bosonic degrees of freedom, second order 
bosonic partial derivatives would appear in eq. (ffTjl and the Schrodinger equation of a 
spinning system of two degrees of freedom would turn out. 

Due to eq. (jMj) . we can make the identification a = 1/h. Hence, taking into account 
fl50|) ■ we have the following relation 

lap 

^ = (72) 

S3 
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